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Abstract. We analyze the time evolution of a one-dimcnsional 
quantum system with zero range potential under time periodic 
^ I parametric perturbation of arbitrary strength and frequency. We 

■ show that the projection of the wave function on the bound state 

I vanishes, i.e. the system gets fully ionized, as time grows indefi- 

Csl ■ nitely. 

o 
o 
o 

I ■ 1. Introduction and results 

^ \ The ionization of atoms subjected to external time dependent per- 

turbations is an issue of central importance in quantum mechanics 
- which has attracted substantial theoretical and experimental interest 

[0) [H- There exists by now a variety of theoretical methods, and a vast 
\ amount of literature, devoted to the subject. Beyond the celebrated 

Fermi's golden rule, approaches include higher order perturbation the- 
ory, semi-classical phase-space analysis, Floquet theory, complex di- 
lation, some exact results for small fields and bounds for large fields 
and numerical integration of the time dependent Schrodinger equation 



0-|jTj]. Nevertheless there is apparently no complete analysis of the 
ionization of any periodically perturbed model with no restrictions on 
the amplitudes and frequencies of the perturbing field. This is not so 
surprising considering the very complex behavior we find in even the 
most elementary of such systems. 



-"^Also Department of Physics. 
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In the present paper we show rigorously the full ionization, in all 
ranges of amphtudes and frequencies, of one of the simplest models 
with spatial structure which, with a different perturbing potential, is 
however frequently used as a model system 0, 0, |jl3[. The unper- 
turbed Hamiltonian we consider is 



d^ 

(1.1) Ho = -—-r^ - gS{x), g > 0, -oo < x < oo. 

Ho has a single bound state Ub{x) = A/poe~^°'^', po = with energy 
—hwo = —h?pl/2m and a continuous uniform spectrum on the positive 
real line, with generalized eigenfunctions 

u(k, x) = ( e^^^" Po^mx\\ _oo < A; < oo 

and energies tx^k"^ /2m. 

Beginning at t = 0, we apply a parametric perturbing potential, i.e. 
for t > we have 



(1.2) H{t)=Ho-gr,{m^) 

and solve the time dependent Schrodinger equation for ifj{x,t), 
^{x,t) = e{t)ub{x)e''^''' 

poo 2 

(1.3) +/ Q{k,t)u{k,x)e-'^^dk (t > 0) 



with initial values 6{0) = 1, Q{k,0) = 0. This gives the survival prob- 
ability |^(t)P, as well as the fraction of ejected electrons \Q{k,t)\'^dk 
with (quasi-) momentum in the interval dk. 

This problem can be reduced to the solution of an integral equation 
ra. Setting 



(1.4) e{t) = l + 2iJ^Y{s)ds 

(1.5) e(A;, t) = 2\k\/[V2^{l - t\k\)] /q Y{s)e'^^+^'>ds 

Y(t) satisfies the integral equation 



:i.6) 



Y{t) = r]{t) |l + ^ [2z + M(t - t')]Y{t')dt''^ 



riit)(l + {2i + M)*Y 
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where h, 2m and | have been set equal to 1 (implying po = 1, uq = 1), 
M{s) = — / — — — 5 — du = — ^ / -^du 

and 

f*9=[ f{s)g{t-s)ds 
Jo 

Theorem 1. Whenri{t) = rsinut the survival probability \0{t)\'^ tends 
to zero as t ^ oo, for all uj > and r 7^ 0. 

Note: For definiteness we assume in the following that r > 0. 
The method of proof relies on the properties of the Laplace transform 
of Y, y{p) = CY{p) = e-P''Y{t')dt' (note that y{p) = ^{l-pCO)). 
In particular we need to show that y{p) is bounded in the closed right 
half of the complex p— plane. Before the proof we describe briefly some 
additional results on this model system, cf. [0. 

1.1. Further results not proven in the present paper. (1) The- 
orem |l| generalizes to the case when ri{t) is a trigonometric polynomial: 

K 

(1.7) r]{t) = sm{jut) + Bj cos{jut)], 

i=i 

where we assume \Ak\ + \Bk\ 7^ 0. 

(2) The detailed behavior of the system as a function of t, u, and r 
is obtained from the singularities of y{p) in the complex p-plane. We 
summarize them for small r; below | < 5 < 1. 

At p = {inuj — i : n ^ Z}, y has square root branch points and y is 
analytic in the right half plane and also in an open neighborhood A/" of 
the imaginary axis with cuts through the branch points. As |Q'(p)| —>■ 
cxD in A/" we have \y{p) \ = 0{ru\p\~^). If |a; — - | > const„0(r^~''), n G 
Z+, then for small r the function y has a unique pole Pm in each of 
the strips —ma; > '^{p) + 1 ± 0(r^^'^) > —niLU — 00, m G Z. ^{pm) is 
strictly independent of m and gives the exponential decay of 6. After 
suitable contour deformation of the inverse Laplace transform, 6 can 
be (uniquely) written in the form 

00 

(1.8) 9{t) = e-^('^'")*F^(t) + e(™"-*)*/i^(t) 

m=— 00 

where is periodic of period 27tuj^^ and 
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Figure 1. Plot of log^o vs. time in units of ujq^ 
for several values of u and r. The main graph was calcu- 
lated from ( |1.8D and the inset used numerical integration 
of (O). 



hmit) ~ ^ Cmjt ^ as t OO, arg(t) G (~2~^'2~^7 
j=o 

A^'o^ too c/ose to resonances, i.e. when Itu — > ©(r^^"^), for all 
integer n, \F^^{t)\ = 1 ± O(r^) and its Fourier coefficients decay faster 
than rl^'"l|m|~l™l/^. Also, the sum in ( pTsD does not exceed OirH'^/"^) 
for large t, and the /i^ decrease with m faster than r'™'. 

(3) By (|1.8| ), for times of order l/F where F = 23^(7), the survival 
probability for oj not close to a resonance decays as exp(— Ft). This is 
illustrated in Figure 1 where it is seen that for r < 1/2 the exponential 
decay holds up to times at which the survival probability is extremely 
small, after which |6'(t)p = 0(t"^) with many oscillations as described 
by ( |1.8| ). Note the slow decay for uj = .8, when ionization requires the 
absorption of two photons. 

(4) When r is larger the polynomial-oscillatory behavior starts sooner. 
Since the amplitude of the late asymptotic terms is O(r^) for small r. 
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increased r yields higher late time survival probability. This phenome- 
non, sometimes referred to as atomic stabilization |jl3|, can be as- 
sociated with the perturbation-induced probability of back-transitions 
to the well. 

(5) Using the continued fraction representation (|2.12| ) F can be cal- 
culated convergently for any uj and r. 

The limiting behavior for small r of the exponent F is described as 
follows. Let n be the integer part oi uj~^ + 1 and assume uj~^ ^ N. 
Then we have, for T > (t = r~'^"'T), 



;i.9) f = -T-Him\n\e{r''-Tf ^"'"^'^^'^ ~ ^ 



' ' ' ' ntu JJ (1 - Vl - muj? 



(6) The behavior of F is different at the resonances ^ G N. For 
instance, whereas if u is not close to 1, the scaling of F implied by ( |1.9|) 



is when cu > 1 and when | < c<j < 1, by taking cu — 1 = j \pl we 
find 

91/4 93/4 

-T-^ lim In ^(r^^T) = 

r^o I ^ '\ 8 16 

a;=l+r2/v^ 

2. Proofs of Theorem |1| 

Lemma 2. (i) CY exists and is analytic in the right half plane M = 

{p : 3fJ(p) > 0}. Furthermore, y{p) ^ as '<s{p) ±oo in H. 

(a) The function y{p) satisfies (and is determined by) the functional 

equation 

(2.1) i/ = r(r- -r+)(/ii + /i2i/ 

with 



(T^f)ip)=f{p±tuj), h,{p) = -^ and h2{p) = ^[l + ^/T^) 

and by the boundary condition y{p) — > as Q{p) — > ±oo in H. 

The branch of the square root is such that for p G H, the real part of 
y/1 — ip is nonnegative and the imaginary part nonpositive. 

Proof, (i) The time evolution of ip is unitary and thus | {ip \ Ub) \ = 
\0{t)\ < 1. The stated analyticity is an immediate consequence of 
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the elementary properties of the Laplace transform 0. The asymptotic 
behavior follows then from the Riemann-Lebesgue lemma, 
(ii) We have in H, 

2„- poo pco 2 

(2.2) £M = lim- / dxe-^'^ / dw 

aio n Jo Jo l + u^ 



^^■^^ 7rl^(l + t.2)(p + ,(l + «2))d« 

For 3?(p) > we push the integration contour through the upper half 
plane. At the two poles in the upper half plane v? + 1 equals and ip 
respectively, so that 



I / (—1) f ds Mq / ds\ i Uq 

n\{2i){p)J s {ip){2iuo) J s) p p 

where Uq is the root of p + z(l + u^) = in the upper half plane. Thus 



(2.5) CM = -'- + 

p p 

with the branch satisfying \/l — ip — 1 as p — in H. As p varies in 
H, 1 — ip belongs to the lower half plane —iM and then a/1 — ip varies 
in the fourth quadrant. 

For 9fJ(p) > 0, u; > we have 

(2.6) C (c^-m) = + 

V / p^ iuj p^ iuj 

(with ^/Y^^ip^^^ = —i^/u — 1 + ip if a; > 1) 
and relation ( ^.1| ) follows. □ 



After the substitution y{p) = 2{y/l — ip — l)e 2^f(p) we get 
(2.7) v{p — iu) + v{p + iu) = -(a/1 — ip — l)f (p) + 



r ' ' ' ' ^- ' u;2 _|_ ^2 



^See also the appendix for a proof of analyticity for 3f?(p) > po (all that is re- 
quired in the subsequent analysis) , relying only on the properties of the convolution 
equation. 
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Remark 3. It is clear that the functional equation ( \2. 1\ ) only links the 
points on one dimensional lattice {p + iZtu}. It is convenient to take 
Po such that p = po + inu with 3?(po) = ^{p) o^nd 

(2.8) 3(po)e[0,a;) 



and write v{p) = v{po + into) = Vn which transforms (\2. iD to a recur- 
rence relation: 



2 / iuj 

(2.9) Vn+i + Vn-1 = -{^/l-^Po + nuJ - l)Vn + 



, up' + (po + inojY' 

where w„ depends parametrically on Pq. It will he seen that the asymp- 
totic conditions as well as analyticity in po determine the solution of 
uniquely. 

Remark 4. The approach is based on a discrete analog of the Wron- 
skian technique. The regularity of the bounded solution of ( ^.dj) will be 
a consequence of the absence of a bounded solution of the homogeneous 
equation 

2 

(2.10) Vn+i + Vn-1 = -(a/1 - tpo + nUJ - l)Vn = DnVn 

r 

a problem which we analyze first. 

Proposition 5. For po satisfying l \2.di ) and 9ft(po) > (actually for 
any po ^ M. = M. U iM.) there is no nonzero solution of i\2.1C\ ) such that 

V G /2(Z). 

Proof. To get a contradiction, assume f ^ is an hC^) solution of 
( |2.1CI| ). Multiplying ( |2.1CI| ) by v^, and summing with respect to n from 



— oo to +00 we get 

oo oo 
(2.11) ^ Vn+lVn+ ^ Vn-lVn 
n=—oo n=—oo 

= 2 ^ ^{vnVn+i) = ^ -y^/l-ipo + nuj -1 



oo 2 



w 



2 



n=— oo 



For G EI the imaginary part of a/1 — ipQ + nu is nonpositive, by 
Lemma ^ and is strictly negative for n < large enough. Thus if for 
some such n, f„ is nonzero then the last sum in (|2.11|) has a strictly 



negative imaginary part, which is impossible since the left side is real. 
If on the other hand f„ is zero when n is large negative, then solving 
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(|2.1CI| ) for Vn+i in terms of the f„, it would follow inductively that 
V = 0, contradicting the assumption. 

□ 

Lemma 6. (i) There is, up to multiplicative constants, a unique pair 
of solutions andv^ of ( 2.10 ) such thatv^ —* as n ^ ±oo ( respec- 



tively). These solutions are related to convergent continued fractions 
representations : 

(2.12) ^± /^;±=:^ = Z}„- 1 



i=Fl/ ■ ± " 1 



Dn±2 

(a) We have the following estimates 



(2.13) -^ = ^ + 0(72-3/2) ^n^±oo) 

where 



Pn Pn 



1 2 , 2 r2-2 + 2ipQ r 

(2.14) — = -V^ (,^>o) 

r r Zry'nuj luon 

1 2i r-— 2 (2-r^)i + 2po r 

Pn ^ r 2r,/\n\uj 2^|^| 

Let be solutions of the one step recurrences = v^^ript- Then 



1 

(2.15) Inf^ = — n\nn + nln I -W — 
2 y 2 y 

+ 2W-+ ( ^" )lnn + o(l) (n ^ oo) 



2ipo + + w\ , I I , . ^ 

+ m \n\ + oil) \n — oo 



(2.16) ln(i;„ ) = — \n\ In \n\ + \n\ In ( ~\ ~ ) + ^7r|r?,| — 2iA/|n|/ci; 
2 \ 2 V f-i^ 

_j_ J. 

and, for some constants , 

(2.17) Hvi)=Hvt) + K^ + o{l] 
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(v^ decay roughly as 1/ \/\n\\ for n ±00, respectively) . 

(in) Two special solutions of ( \2.1C\ ), and v~ , are well defined by: 



(2.18) 



Y[ -q: for n > N, and f „ = Y[ ~~ f^f^ ^ < ~^ 

j>n+l Pi j<n-l 



if N is sufficiently large (this amounts to making a convenient choice of 
the free multiplicative constant in (i)). These functions do not depend 
on N . and v~ are linearly independent for po G H; their discrete 
Wronskian, defined by W{v~^ ,v~)n = v^v~_^_i — v~v^_^i, satisfies 



(2.19) W{v+,v-) = const j^O 

As functions of parameters, and W{v~^ ,v~) are analytic in Pq G H. 
If oj ^ {0,72^^ : n G N} then and W{y^,v~) are analytic in some 
neighborhood of po = as well. For any u > 0, are Lipschitz 
continuous of exponent at least 1/2 in po, forpo G R. 

Proof, (i) We look at v~^, the case of v~ being similar. Dropping the 
superscript we have from ( |2.1U| ) 



(2.20) p„ = 

J^n — Pn+1 

To find the analytic properties of the solution it is convenient to 
regard ( |2.2(J| ) as a contractive equation in the space i°°{SN) of sequences 



{pj}j>N in the norm ||p||oo = sup^-^^ \pj\. Let be large. The map 
J : Sn ^ Sn defined by 



(2.21) J(p), 



Dn ~ Pn+1 

depends analytically onpo and is Lipschitz continuous of exponent 
at least 1/2 if 3?(po) ^ 0. In addition, if ||pj||oo < 1 we have for 
sufficiently large = N{pq, lj, r) 

(2.22) II J(p)lU < . ^ ^ = ^ < 

f [^/\Nuj\ - 1 - bol - 1) - 1 

Similarly, 





r 


1 


\U! 
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(2.23) 

\\J{P)-J{P') 



< 



IP-P lloo 



\po\ - 1 



< 



N\uj\ 



\P- P lloo 



for sufficiently large which shows that J is contractive in the unit 
ball in I'^^Sn)- Thus, equation ( p.21|) has a unique solution in S'at, 
which depends analytically on po ^ H and is Lipschitz continuous of 
exponent at least 1/2 if 3?(po) > 0. This also implies the convergence 

of (iH). 

Note that given /Ci C H and K,2 C both compact, can be 
chosen the same for all po ^ and r G /C2- 

(ii) From ( ^.22| ) it is seen that \pj\ = 0{j^^^^) for large j. Thus, we 
may write, for large j, 



(2.24) 4 = 



The estimates ( |2.13|) now follow by a straightforward calculation. Since 
In f + = \'nv~^ + X]j'=iv+i In p/) the estimates follow from ( p.24|) and the 
Euler-Maclaurin summation formula. 

(iii) As before, we only need to look at v~^. We take two compact sets 
/Ci and /C2, and choose A^ as in the note at the end of the proof of (i). 
Taking the log in the definition ( p.l8| ), the infinite sums are absolutely 
convergent. By standard measure theory, f + has the same analyticity 
properties in the interior of /Ci x /C2 and Lipschitz continuity in /Ci x /C2 
as those of p"*", when n > N. Now, ( p.lO| ) easily implies that the same 
is true for < A^ as well. 

If fn and gn are solutions of (|2.10|) then {gn+i + gn-i)fn - (/n+i + 

fn~l)gn = and thus Wn{f,g) = fngn+l - gnfn+1 = COUSt. ThuS, if 

W^ra(/? fi') = foi' some n then Wn = and / = const g. The smoothness 
properties follow from the proof of (iii). 

□ 



Proposition 7. There exists a unique solution of Tj ) which is bounded 
as Q'(p) ±cxD in H. This solution is analytic in p E M, and v{p) = 
0{p-^) as Q{p) ±00, p e H. 

Proof. By analyticity and continuity W{v^, v^) does not vanish for any 
p G EI and r > 0, a; > 0. By Lemma |Tl| the function v defined through 

v{po + inuj) = fn, where 
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(2.25) f,:=W{v^,v-r' U+ Yl vfH, + v-J^^^^, 

\ l=—oo l=n 

and 

(2.26) H„ ^ ^^^^ 

has the required properties. Since no solution of the homogeneous 
equation is bounded on Z, v is the unique solution with the desired 
properties. □ 

Note 8. The link between y and fn is 



(2.27) y{p) = 2{^l-tp-l)e-2^U; for p = po + 



znuj 



Proposition 9. The function y{p) is analytic in the right half plane, 
Lipschitz continuous of exponent at least 1/2 on the imaginary axis and 
limp^o y{p) = i/2. 

Proof. Since W is analytic in H, continuous and nonzero in H, W is 
bounded below in compact sets in H. Then, the smoothness proper- 
ties of y derive easily from those of := W{v~^ ,v~)fn on which we 
concentrate now. 



(a) For n >2 we write, using ( p.26|) . 



re-l 



(2.28) qn = v^Yl ^THi + ^ 



/=— cxD l=n 



+ v^iue 2^ ; — — + 



,Po(po-2iu;) po{pQ + 2iuj) 
The last term in parenthesis can be rewritten, using also ( p.lO| ), as 



, iiv^ — f i) 2uj /"Wi + f 1 

(2-29) \i , + ' 



pI + 4cj2 pI + 4^0-2 \^ 



Auj ^/l - ipo - 1 _ 



pI + 4cj2 r{pl + Aoj'^) Po 



12 



O. COSTIN, J. L. LEBOWITZ AND A. ROKHLENKO 



Thus we see that g„ is continuous as 3ft(po) ~^ and Q{po) G [0,ti;) [cf. 



(|2.8| )], if n > 2. A very similar calculation shows the continuity of g„ 
if n < -1. 

(b) By part (a), y{p) is continuous as 3?(p) J, with '^{p) > 2 or 
Q'(p) < 0. Now, (|2.1|) written in the form 



(2.30) rph2{p)y{p) = rp[ hi{p + 2iuj) — hi{p) 

_l_ rp(y{p + iuj) + h2{p + 2iuj)y{p + 2iuj) 



shows that y{p) is Lipschitz continuous as 3?(p) i if ^{p) > —2 thus 
for all Q{p). The value of y{0) is easily calculated using ( |2.30| ). □ 

Proposition 10. I + 2i lim^^^oo Jq Y{s)ds = 0. 
Proof. Indeed, 

poo 

(2.31) 27rz / Y{s)ds 



'0 



lim lim ( / + / — (i/2 + (y(p) - i/2))dp = 

x^ooS^O+ VJ-ioo JiS ) V 

□ 

3. Appendix 



Lemma 11. Equation (lj_6) has a unique solution Y G Ll^^{M.'^) and 



\Y{x)\ < Ke'^'' for some G M+ and C G 

Proof. Consider L]^c\^, A\ endowed with the norm := \F{s)\e~'^'^ds, 

where > 0. If / is continuous and F,G E LloA'^^ ^ straightforward 
calculation shows that 



(3.1) ||/F||.< ||F||,sup|/| 

[0,A] 

(3.2) 11^*6-11, < 

(3.3) \\F\\u —^0 as u ^ oo 

where the last relation follows from the Riemann-Lebesgue lemma. 
The integral equation ( |1.6|) can be written as 

(3.4) Y = rr] + JY where JF := rr]{2i + M) * F 
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Since M is locally in V" and bounded for large x it is clear that for 
large enough C2, and for any A, (|1.6| ) is contractive if z/ > C2. □ 
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